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Gravitational waves generated by the final merger of compact binary systems depend on tfie 
structure of tlie binary's members. If tfie binary contains neutron stars, measuring such waves can 
teach us about the properties of matter at extreme densities. Unfortunately, these waves are typically 
at high frequency where the sensitivity of broad-band detectors is not good. Learning about dense 
matter from these waves will require networks of broad-band detectors combined with narrow-band 
detectors that have good sensitivity at high frequencies. This paper presents an algorithm by which 
a network can be "tuned", in accordance with the best available information, in order to most 
effectively measure merger waves. The algorithm is presented in the context of a toy model that 
captures the qualitative features of narrow-band detectors and of certain binary neutron star merger 
wave models. By using what is learned from a sequence of merger measurements, the network can 
be gradually tuned in order to accurately measure the waves. The number of measurements needed 
to reach this stage depends upon the waves' signal strength, the number of narrow-band detectors 
available for the measurement, and the detailed characteristics of the waves that carry the merger 
information. Future studies will go beyond this toy model, encompassing a more realistic description 
of both the detectors and the gravitational waves. 

PACS numbers: 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 



Much of the promise of gravitational-wave (GW) observation is in its potential as a novel probe of physics and 
astrophysics. Because GWs couple very weakly to matter and arise solely from gravitational interactions, sources that 
are completely dark electromagnetically may be strong GW emitters. By tracking and measuring the waves generated 
in violent astrophysical events, we may gain insight into physical processes that cannot be easily measured in other 
ways. 

One example of such hard-to-measure physics is the late merger of compact binary systems. GW emission carries 
energy and angular momentum out of the binary, driving the compact bodies ever closer together. During the early 
inspiral portion of this coalescence process, the structure of the binary's members plays little role; they can be usefully 
approximated as point masses, or spinning point masses. As the bodies come closer together, their internal structure 
becomes very important. The GWs generated in the final stages of inspiral and merger, when the bodies collide and 
merge into some coalesced state, will carry information about the bodies' structure. 

It has long been recognized that, if the binary contains at least one neutron star, merger waves will depend on the 
nature of neutron star matter (cf. Refs. |l| 3, ^ and references therein). This opens the exciting possibility that 
GW measurements by detectors such as LIGO [5] could study the properties of very dense matter, such as its equation 
of state (EOS) |[ 0|. Of particular interest will be testing whether neutron stars contain a core of "exotic" matter, 
such as a free quark state of some kind. Models of compact stars comprised of a free quark fluid, or with a free quark 
core, can have a structure rather different from ordinary neutron stars — they are often of smaller radius, and there 
may exist a sharp density transition at some finite radius from the star's core ^, ^- If one or both members of a 
binary system had such a structure, there could be an observable imprint on the merger GW signature. (We should 
emphasize, though, that much work remains to evaluate whether such a signature exists, and if so, what is its nature.) 

Simulations of binary neutron star merger using Newtonian gravity and a polytropic EOS (see, e.g., jl^, [ll)) have 
found that the merger waves indeed carry information about the EOS, but that these waves are at very high frequency 
(/ ~ 1500 — 3000 Hz) where broad-band LIGO-type detectors do not have good sensitivity. More recent work using 
irrotational matter configurations in the conformal approximation to GR shows that the EOS-dependent structure 
is likely to come out at somewhat lower frequencies, / ~ 1000 Hz; this is still high enough that broad-band detector 
sensitivity is not very good. If one member of the binary is a black hole, the EOS-dependent information will come 
out at still lower frequencies, / ~ 400 — 1000 Hz — the larger system mass shifts all frequencies downward. In this 
case, broad-band detectors are more useful for studying the merger waves, but still may not be ideal. 
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To learn as much as possible from the waves generated during the merger, broad-band GW detectors should be 
supplemented by narrow-band detectors. Acoustic detectors in existence today |l^ and special interferometer 
topologies under development such as signal recycling and resonant sideband extraction have good sensitivity 
in a narrow band at high frequencies, 5/// ~ 0.1 — 0.2 for / ~ 500—2500 Hz. Narrow-band detectors answer essentially 
a yes-no question: "Did the binary radiate in my frequency band?" A "xylophone" of narrow-banded detectors would 
probe gross features of the merger waveform, such as a sharp cutoff in the GW spectrum (seen in recent simulations 
1 12 1 ) , or the waves generated by a transient bar that forms in the merger detritus (seen in some Newtonian simulations 
1 10, Such simple measurements should be robust in the sense that they wouldn't require detailed modeling of the 
waves' phasing — very important, since theoretical uncertainites in the merger waveforms are likely to be significant 
even when these measurements can be made. 

Practical considerations such as cost and available facility space will limit the number of narrow-band detectors 
that can be used for each measurement. To make best use of these detectors, the network should be designed in a 
way that is in some sense optimal: the narrow-band detectors should be configured, in concert with the broad-band 
detectors, so that the network of all detectors is most likely to provide new information about merger waves, given our 
best present knowledge of the waves' properties. How one "tunes" a detector network in this manner is the subject 
of this paper. 

We assume that an inspiral has already been measured, so that we know merger waves must be present in the data. 
For the purposes of this analysis, we assume further that the waveform depends on a single parameter A that grossly 
characterizes the merger waves. This A could be the frequency of a sharp cutoff in the wave spectrum, or the frequency 
at which a short-lived bar may radiate for several cycles. Theoretical modeling allows us to phenomenologically relate 
this parameter to a description of the binary's stars. For example, in models in which the wave spectrum sharply cuts 
off, A is most strongly related to the compactness of the stars: smaller stars exhibit a cutoff at higher frequencies 
since they spiral in further before the spectrum cuts off. Describing the merger features with a single parameter is no 
doubt an oversimplification, but is useful for demonstrating how network tuning works to zoom-in on gross features 
of the waves. Tuning the network means finding the configuration which measures A with as little error as possible. 

We develop a tuning algorithm that does just this: it configures the network to measure A as accurately as possible, 
given our uncertainty in A's value, and updates ("retunes") the network as measurements teach us about merger 
waves. This algorithm is based on the maximum likelihood GW measurement formalism developed by Finn [g_7|. 
Finn defines two probability distributions which play a major role here: the prior probability, po(A), summarizing all 
that is known about A before measurement; and the posterior probability, Ppost(A|A), summarizing what is known 
afterward. The posterior distribution is built from the measured datastream, and so explicitly depends upon the 
detector network's characteristics and on the measured waveform h{X), where A is the unknown, true value of A. 

Following measurement, the posterior probability is the tool one uses to estimate the values of the parameters 
describing GWs, and also to estimate the error in those values ^ It is not quite the tool needed here: we 
need to estimate the accuracy with which A is likely to measured, but we need this estimate before measurement. 
To this end, we introduce an additional probability, the anticipated distribution of A. This is a marginal distribution 
found by integrating out the dependence of the posterior probability upon the unknown true parameterization A: 

Pa„t(A)= J po{X)Ppost{X\X)dX. (1.1) 

This is the distribution that we we anticipate will describe A after a measurement. Like the posterior probability, it 
depends on the detector network, so we can tune the network's adjustable parameters to find the network which we 
anticipate will measure A as accurately as possible, given our current ignorance of the merger waves. 

From the anticipated distribution's definition, it is simple to update and improve the network as we learn more 
about A. Following a measurement, we construct the posterior distribution Ppost(A|A) from the datastream. We then 
use this posterior distribution as the prior distribution for the next measurement: Po(A) — Pp~^(A|A). We update 
the anticipated distribution with the new priors, and then update the network. The detector network is thereby 
adjusted and improved following each measurement, so that our ability to measure A is incrementally improved 
by each measurement. Our gradually improving knowledge of A's value is manifested as a gradual peaking of this 
distribution: we begin with a prior describing complete ignorance (uniform distribution) and find that after several 
merger measurements the distribution begins to peak around the true value A. The more we learn about the merger 
waves, the more sharply peaked becomes this distribution. 

In this paper, we demonstrate the concept and principles of network tuning using a toy model for the merger and for 
the narrow-band detectors. In Sec. we describe in more detail the probability distributions introduced above, and 
then discuss how we use them to tune our network designs. Our tuning procedure is given explicitly at the end of this 



section. We next present our toy model in Sec. [H: we approximate narrow-band detectors as zero bandwidth (delta 



function) GW detectors with adjustable center frequency, and treat the waveform h{t; A) as a quadrupole inspiral 
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chirp up to a merger frequency fm = A x 1000 Hz. This highly simpUfied description throws away important physics, 
particularly the finite bandwidth expected in real narrow-band detectors and the slower frequency roUoff seen in recent 
merger computations (e.g., Ref. (l^). It is simple enough, though, that many of the calculations needed can be done 
analytically, and is close enough to the real problem that it should give a good sense of how network tuning is likely 
to proceed in practice. 



We test our tuning algorithm in Sec. IV. We show that measurements can converge on an accurate value for A 
after measuring some number of binary merger events. The size of that number depends on the signal strength (a 
few strong signals can drive convergence rather quickly) and the number of narrow-band detectors available for the 
measurement (having at least two available can speed up convergence quite a bit). The good behavior of our tuning 
procedure can be taken as an indication that this algorithm is robust. However, we cannot pretend that our analysis 
is in any way definitive: the toy description of the detectors and the waveform neglects several important effects. This 
analysis should therefore be regarded as a proof-of-concept presentation; future work, discussed in Sec. will put the 
various complications neglected here back where they belong. In particular, we plan to use more realistic descriptions 
of the narrow-band detectors (cf. Refs. ^), and merger waveforms taken from recent computational models. 

II. FORMALISM 
A. Probability distributions 

As mentioned above, we assume that we have already detected a GW via the inspiral, and so we have already 
learned many interesting parameters characterizing the binary — combinations of the masses and spins, the distance 
to the binary to some accuracy, etc. We now wish to measure the merger parameter A. 

Consider a network of N GW detectors. We use the notation of Appendix A of Ref. : the datastream of the 
j-th detector in the network is gj] the output of the whole network is g. Because we assume that a GW has been 
detected, this datastream consists of a GW signal h{X) plus noise n: g — h{X) + n. (Recall that A represents the true, 
unknown value of A that describes the waves in the datastream.) Our goal is to determine the probability that, given 
this datastream, we will measure the parameter A. 

Our first step is to find the posterior probability that we measure the waveform h{X) given the datastream g. Finn 
has developed this probability for a single detector ; using the notation of Ref. , the generalization to multiple 
detectors is straightforward. Following Finn, we begin with Bayes' law: 

P[MA)|,1^M»M. (2.1) 
P{9) 

The probabilities written here have the same meaning as those used in Eq. (2.3) of JT?! : 

P[ft,(A)|(7] = the conditional probability that h{X) is present given the datastream g ; 
P[(7|/i(A)] ~ the conditional probability of measuring the datastream g given h(X) ; 
P[h{X)] = the prior probability that the waveform h{X) is present 
= Po(A) ; 

P{g) = the probability that the datastream g is observed 

^ 1 . (2.2) 

We have simplified the last two of these distributions. The first, P[h{X)] = po{X), follows because by assumption 
the waveform depends only on A. This is no doubt an oversimplification — other parameters, such as the mass ratio 
of the binary and its spins, are likely to have an important influence as well. To keep this analysis simple, we will 
focus on a single parameter model for the merger; this is adequate for now, and perfectly formulated for the toy 



model presented in Sec. HI. The second simplification, P{g) = 1, follows because g has already been measured and is 
therefore determined. Here we depart a bit from the analysis of Ref. jT^ — Finn does not assume that this particular 
datastream g has already been measured. His g is undetermined, and has a non-trivial probability distribution. 

The distribution P[g'|/i(A)] used here is identical to that used by Finn (modulo the extension to multiple detectors). 
It answers the question "What is the probability that, if the data contains the GW h{X), noise is present such that 
n + h{X) is consistent with the observed data stream g?" This shows that P[g\h{X)] is equivalent to P[g — ft.(A)|0], 



the probability that g — h{X) is pure noise. Under the assumption that the noise is Gaussian, this is p7| 



P[g\h{X)]^ICexp 



-Kg-km-hw) 



(2.3) 
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~aAfr[SkifrYhif)df , 



where the inner product is given by 

(d\b) = 4Re 
^ ^ JO 

and /C is a normalization constant. Here, aj{f) is the Fourier transform of aj{t): 



(2.4) 



(2.5) 



and the asterisk denotes complex conjugation. The matrix Sh{f) is a generalization of the (one-sided) strain noise 
spectral density. Diagonal components S^(/)jj represent the usual spectral density for detector j; off-diagonal elements 
represent correlations between detectors j and k. For more details, see Appendix A of Ref. |Q. Combining these 
distributions, using g — h(X) + n, and defining Sh = h(X) — h{X) yields 



P[h{X)\g]=ICpo{X)exp 



l{6h\Sh)~{n\6h)~^{n\n) 



(2.6) 



Equation ( p.6[ ) depends on the specific noise instance n{t) pre sent in the detector network. We circumvent this by 
taking a "frcqucntist" viewpoint and ensemble averaging Eq. ( |2.6| ) with respect to the noise. The resulting distribution 
is our posterior distribution for A, given data containing A: 



Pr 



post 



(A|A) ^ po{X)E[p[hiX)\g]}_^ 

= poiX) f P{n)P[h{X)\g]Vn 



ICpo{X) / Vn exp 



-- (r?|n) 



exp 



[6h\Sh) - {n\Sh) - 1 in\n) 



(2.7) 



The final integral is a functional integral, taken over all possible network noise instances n. We evaluate it in Appendix 
^; the result is 



-Ppost(A|A) = /Cpo(A)exp 



--iSh\Sh) 



(2.8) 



This distribution is normalized so that J Ppost(A| A)(iA = 1. The normalization constant /C thus depends on A and the 
noise characteristics. 

As discussed in the Introduction, i-'post(A|A) is not useful for designing a detector network since it is constructed after 
measurement. This is reflected by its implicit dependence on A. A more useful quantity is obtained by marginalizing 
on A. This gives a distribution which represents the posterior distribution that we anticipate will describe the A after 
measurement, modulo our ignorance about A: 



Pant (A) = / po(A)Ppost(A|A)dA. 



(2.9) 

The anticipated distribution answers the question "How is A likely to be distributed given our ignorance about the 
neutron star merger waveform?" 



B. Tuning the detector network 



We now choose the detector configuration that we anticipate will measure A as accurately as possible. Consider the 
mean and variance of the distribution ( ^.9| ): 

Aant = J A Pant (A) dA, (2.10) 
'^ant = / (A- Aant)'Pant(A)dA, 
A^ Pant (A) dA — Aanf 



(2.11) 
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The m ean is indicative of the typical value of A that a particular observation will measure. Likewise, the variance 



11) is an indicator of how accurately A will be measured; we anticipate that a measurement will typically find A 



ant ■ 



within (Tant of A 

Through the inner product (2^), the variance ( 2.11 ) depends on the characteristics of the detector network. Some 
of these characteristics will depend on adjustable parameters, such as the central frequencies of narrow-band detectors. 
Let the kth such adjustable parameter be fj,k. The network which minimizes dant is the one which we anticipate will 
measure A most accurately: 

miuCTant The optimally tuned network . (2-12) 

We call the process of finding the parameters which minimize (Tant "tuning" ; the network so produced is "optimally 
tuned" or simply "optimal" . 

Because of the non-Gaussian nature of the relevant distributions, a variance such as cra,jt is a less than ideal measure 
of error. If we were truly interested in characterizing the error of a measurement, we would instead study contours of 
a likehhood function (see, e.g., Refs. 0, |l^). However, a^^^ is very easy to compute, and has the correct qualitative 
behavior — it will get smaller as more rigorous measures of error decrease. Tuning our network in such a way that 
CTant is minimized will drive more rigorous notions of error to be small as well. 

Since Ppost(A|A) is built from the measured data, it depends on the measurement's signal-to-noise ratio (SNR). 
Thus the anticipated probability, and the optimal network it produces, likewise depend on SNR. This is arguably a 
weakness of this prescription for network tuning — the network found by minimizing a'^^^ depends explicitly upon the 
anticipated strength of the waves that we hope to measure. For our present purposes, we believe that leaving the SNR 
in the algorithm as an adjustable parameter is actually desirable. Properly eliminating the SNR dependence would 
require integrating it out with a correctly constructed probability function for its distribution. Although there exist 
estimates of the compact binary coalescence rate which we could use to construct such a function, these estimates 
are currently rather uncertain By the time that the interferometer configurations we consider here can be 

installed, LIGO-type detectors could be detecting many events per year, or just a few. If many are detected, then 
observers may choose to focus their merger analyses on the strongest events in this large set. On the other hand, if 
the rate turns out to be low, observers will be forced to work with relatively weak signals. It is therefore useful to see 
how SNR choice affects network design, exploring how tuning will work as a function of the SNR levels that turn out 
to be appropriate. Accordingly, we choose to leave in the SNR as a parameter for now. 



Finally, note that the prior probability po(A) plays a crucial role, through Eqs. (2^) and ( |2.9| ), in determining the 
optimal detector network. This gives a natural way to update the network as our knowledge of merger waves improves: 
we continually update the priors as we learn more about A, and concomitantly update the network to zoom-in on the 
true value of A. This updating is the key element of the network tuning algorithm. 

Explicitly, this algorithm works as follows: 

1. Before the first merger is measured, our understanding of A is likely to be crude. Take the prior distribution to 
be uniform between some upper and lower limits: po{X) = const, Amin < A < Amax- 



2. Configure the detector network by minimizing Eq. ( ^.11 ) with this uniform prior 



3. Measure a merger. Build the posterior distribution Ppost(A|A) [Eq. ( |2.8| )] from the measured data. 

4. Use this posterior distribution as the prior for the next measurement: set Pq''*(A) = -Ppost(A|A). 



5. Reconfigure the network by minimizing (2.11) with this prior. 

6. Go to 3. 



As we will show in Sec. [V, this algorithm effectively zooms in on the true parameterization A of the waves in the 



datastream, at least for the toy model we present next. 



III. TOY MODEL 



Here we describe the simplified model we use for the waveform (Sec. HI A ) and for t he de tector noise (Sec. Ill B| ) . 
These models make it possible to eval uate t he posterior distribution analytically (Sec. Ill C| ) . We discuss the short- 
comings of this toy description in Sec. HID. We plan to address these shortcomings in future work. 
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A. Waveform 



The gravitational waveform has two polarizations, /i+ and hx (named for the axes of force lines associated with 
them). Each detector in the network measures some linear combination of these polarizations: 

/i, =F,-+/i++F,-x/ix. (3.1) 

The fmictions -Fj,+,x depend on the type of GW detector used for the measurement. For interferometers, standard 
formulae for -F+,x are given in, e.g., Ref. (2| [Eqs. (104a, b) and preceding text]. For certain spherical acoustic detectors 
(cf. Ref. ||l^ and references therein), F+ ~ Fx ~ const — spherical antennae have nearly equal sensitivity to both 
polarizations, independent of the source's sky position. 

The polarizations of the leading quadrupole harmonic of the inspiral waveform are, in the frequency domain, 

h+{f) = 2{l + cos^i)hQ{f), 

hQif) ^ V4 ^ rv^.mf), (3.2) 

V 96 r 

The angle l is between the line of sight to the binary and its orbital angular momentum, r is the distance to the 
source, M. = (mim2)^/^/(TOi + 7712)^^^ is the binary's "chirp mass", and $(/) is a phase function, described further 
below. 

To make our toy model for the waveform, we first simplify this inspiral in two ways. First, we restrict ourselves to 
equal mass binaries: we set mi = m2 = 1.4 Mq. Since we are focusing on neutron star binaries, this is not a drastic 
simplification. Second, we use only the Newtonian, quadrupole contribution to the phase, ignoring higher-order 
post-Newtonian contributions: 

$(/) = 2^fU '<t^c-\ + \ {SnMfr'/' , (3.3) 

where is the time at which the frequency (formally) diverges to infinity, and 0c is the phase at that time. 

To represent the transition from inspiral to merger, we multiply h+^x by 6(A/kHz — /), where 9 is the step function 
and /kHz = 1000 Hz. We assume that A lies between 0.5 and 1.5; this is consistent with the results seen in Ref. [^ . 
This extremely simple waveform captures the essence of coalescence models in which the GWs shut off very rapidly as 
the neutron stars come into contact with one another. Because we neglect higher order contributions to the phasing, 
and because our model shuts off the wave emission far more rapidly than is seen in any actual calculation, our results 
should be considered illustrative rather than definitive. 

The frequency domain GW signal measured in detector j of the network is therefore 

h,{f) = [2(1 + cost2)F,-+ +4coste'"/2F,- x] hQ{f)Q{\hn. - f) ■ (3.4) 
A very useful quantity which we can calculate from this waveform is the SNR pji 



Sh{f)j 



- ^1 ^^rf/ (3-5) 



5|*||7r-4/3A^5/3 |-AAh. f-7/3 



-df . (3.6) 



24 r2 J, Shif) 

The first line is a standard formula for the SNR of a GW measurement (cf. Refs. [0 |l8[ |l|l). The quantity Shif)j 
is the spectral density of noise in the j-th detector; we discuss our toy description of this quantity in the following 
section. In the second line, we have specialized to the toy waveform. We have introduced here the angular function 

\^\^. ^ A [{1 + €03"^ if Fl^+ A cos'^ iFfx] ■ (3.7) 

Note that |*|^ is identical to the 9^ defined in Ref. |l|l; we use a slightly different notation to avoid confusion with 
the step function. This function lies in the range < < 16; its value averaged over all sky positions and source 
orientation angles is = 2.56. The distribution of is quite asymmetric, with values smaller than \'i>\'j more 
likely than larger ones. See Ref. iQ for further discussion. 
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We plot a representation of this waveform in Fig. ^ comparing it to our model of the detector noise (discussed in 
detail in the following section). The waveform plotted here is somewhat massaged — we actually plot 



^iot(/) = ^^f\h{f)\ 



(3.8) 



This quantity has the same units as the noise spectrum, and has the nice feature that the integral of ^piot(/)^ over 
the noise is the squared signal-to- noise ratio p^: using Eq. (|3.5|), we see that 



d\nf 



/iplot(/)j 

Sh{f)j 



(3.9) 



The waveform we show corresponds to a pair of \AMq neutron stars with optimal orientation (|^'P — 16) at 500 
Mpc; or, with "average orientation" (l^'P = 2.56) at 200 Mpc. We have chosen A = 0.8, so that the waveform's 
merger cutoff is at 800 Hz. 



Inserting the toy waveform (p.4|) into Eq. (2.8) yields 



^'post(A|A) = po(A)cxp 



-Re / [2(1 + cosi2)F_,-+ -f'4cosie"/2^j;^][2(i + cosi2)Ffc.+ 
Jo ' 



2 _ 



4cosie^^/2j^fe,x] e (A/kHz -/)- e (A/kHz - /) " |/iq(/)|2 [S„(/)- ^/ 



(3.10) 



Define Alo = min(A, A), Ahi = max(A, A). Because of the functions in Eq. (3.10), the domain of integration is 



restricted to the range /lo < / < /hi (where /lo = Alo /kHz, /hi = Ahi /kHz)- For the range of A that we have 
s elect ed, the detector noise in the merger band is dominated by laser shot noise. We will use this to simplify Eq. 
(3.10) further in what follows. 



B. Detector noise 



Very detailed parameterized models for the LIGO noise budget exist and can be manipulated using codes such as 
BENCH [ p6| (which was used for the analysis of [^l| ). For our purposes, much simpler descriptions that are simple to 
manipulate analytically will be mostly adequate. 

We will take our broad-band detectors to have the noise curve anticipated for the second-generation LIGO detectors 
("LIGO-II noise"), as described in Refs. |2^, In particular, since we focus on measurements of binary neutron star 
inspiral, we will assume the broad-band measurements are made using the interferometer configuration optimized for 
that source (the curve labeled "NS-NS" in Figure 1 of Ref. [Q); we show it in Fig. |l]. 

The NS-NS LIGO-II noise profile we use is a broad-band signal-recycled configuration. In the shot noise dominated 
band, its noise spectrum grows with an /^ law: 



fl.5 X 10"^^ Hz" 



(3.11) 



The numerical value of the parameter 7ns-ns given here was found by fitting to the high-frequency end of the NS-NS 
LIGO-II noise curve given in p^ . 

We will take our narrow-band detectors to be LIGO-II detectors operated in a narrow-band configuration; cf. 
discussion in Refs. [|2^, |2l|. A detailed description of the noise in this configuration can be found using bench; near 
the most sensitive frequency of such a configuration, it is well described by the approximate formula M, E^, Pq| 



Srif) =^ 25o ( ^ 
Jo 



The parameters 5*0 and /o appearing here are 



A/ 



(3.12) 



So = 
fo = 



(3.13) 
(3.14) 
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The numerical values of the laser wavelength Al, laser power /q, photodiode efRciency 77, mirror power loss (3, and 
armlength L are given in Table II of Ref. . The frequency fn is controlled by the position of one of the mirrors in 
the narrow-banded interferometer, and as such is not too difficult to adjust; it might even be possible to dynamically 
adjust this position, in order to track the evolution of a particularly interesting feature. 

The bandwidth A/ is rather more complicated. It is largely set by losses in the mirrors, but also depends in an 
important way on the position of the mirror that controls /_r. Hence, as we tune /ij, the bandwidth changes as 
well. For the purpose of our toy model, we will ignore this effect; examples of the toy narrow-band noise spectrum 
(with fixed bandwidth) are shown in Fig. ^ We emphasize that the bandwidth cannot be treated as constant in a 
realistic analysis. Indeed, there are many coupled effects that appear when a realistical description of a narrow-band 
interferometer is used. For example, configurations that improve detector response in a narrow-band at high frequency 
can actually improve response at low frequency as well p9| , pO| . 

For the toy model, we will approximate the narrow-band detectors as having infinitesimal bandwidth: 

^ -■^5{f-fR)- (3.15) 



This form is particularly easy to use, and is not a bad description of the interferometer response around the center 
frequency when the bandwidth is very narrow. The prefactor is chosen so that the integral of 1/S'^^(/) is the same for 



Eqs. (3.12) and (3.15). We will treat the frequency fu of each narrow-band detector in our network as an adjustable 
parameter, and tune our network by varying those frequencies to find the optimal network for measuring merger 
waves. Note that by using such a narrow detector response, we throw away potentially useful information. This is an 
obvious point of improvement for an analysis that goes beyond this toy model. 



C. Evaluation of the posterior probability 



We now use the toy description of the noise to evaluate analytically the posterior probability distribution (3.10). 
First, we assume that the various individual detectors on the network have uncorrelated noises: we set the off-diagonal 
terms in the matrix Sh{f)~^ to zero. This should be a reasonable assumption in the high frequency, shot-noise- 
dominated regime appropriate to this analysis — each interferometer will have its own laser and optics system, and 
should therefore be reasonably isolated from al l oth er detectors (though events affecting interferometers at a common 
site could violate this assumption). Equation (3.10) then simplifies to 



-Ppost(A|A) = _po(A) exp 



N /./hi 

j — l -'JLO 



■df 



(3.16) 



where N is the total number of detectors, broad-band plus narrow-band. 

To simplify this analysis further, we will assume that all detectors are housed within a single facility, so that the 
functions take a single value jvE'p for all j. Further, we will assume there is only one broad-band detector in this 
facility, and Nn narrow-band detectors. The sum in the exponent of Eq. ( |3.16 ) can now be simplified: 



N 



2 l^""' 



/l 



Shif), 



1*1 



/l 



df 



Afjv „/h 

._i J f 



' \hQi.f)\' 



df 



(3.17) 



Using Eqs. (3.11) and (3.15), these integrals are simple to evalute. The result is 

o2 



-Ppost(A|A) = po(A)exp 



3 P] 



insp 



1 



^-10/3 

/kHz 



fo 



20 ar 7NS-NS 

Nn / j: \ 7/3 



1 



1 



16 (77 



,-7/3 V- 
JkHz / . 



/kRg 

fR,j 



XlO/3 

(A/kHz — fR,j) — (A/kHz 



fR,j) 



(3.18) 



The frequency fuj is the adjustable frequency of peak sensitivity for the j-th narrow-band detector. The step 
functions enforce the fact that the integral over the narrow-band noise is non-zero only if /lo < fR,j < /hi- We have 
simplified this expression considerably by normalizing the amplitude of the GWs to the inspiral SNR Pmsp measured 
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in the broad-band detector. Using Eq. (3.6), 



2 5|^|2 ^-4/3_V(5/3 rcc. j-7/3 

5|^|2 ^-4/3^5/3 



24 



cry . (3.19) 



Taking the upper hmit of the integral to infinity (rather than A x 1000 Hz) is fine for the inspiral SNR since most 
of the signal accumulates near / ~ 100 Hz — the upper frequency has little effect. With the broad-band LIGO-H 
noisecurve used here, cry ~ 8.3 x 10^"* Hz^^^^. For a binary consisting of a pair of 1.4 M© neutron stars, the SNR 
becomes 

M^V'^ /200Mpc\ 



With this parameterization, we see how the tuning algorithm will depend upon SNR, as described in Sec. 11 B . Note 
that, in this toy model, the inspiral SNR and the frequency completely determine the SNR measured in the j-th 
narrow-band detector: 



S^^if), 



Phisp /o .-7/3 



16 CT7 Sa-'^'^ 



Ir ■ (3-21) 



Plugging in the values of /o, 5*0, and cry, we find 



^=0.05fi^V'\ (3.22) 

Pinsp V JR,j J 

The narrow-band SNR drops off very quickly with frequency. This follows directly from our toy model for the waveform 
— the amplitude decays with frequency as /"^^^ up to the step function break, so the SNR decays likewise. 



The anticipated probability distribution is found by plugging ( 3.1g ) into (p.g|) . It is best evaluated numerically. 



D. Shortcomings of this toy description 

Although useful for demonstrating how network tunin g wo rks, we cannot pretend that the toy model is anything 
more than illustrative. First, the waveform we use [Eq. ( |3.2D ] is drastically oversimplified. Although it produces an 
energy spectrum qualitatively similar to that seen in recent computations (e.g., Ref. |12|| ), it shuts off the inspiral 
waves far more quickly than those calculations predict. In the spectra developed in Ref. |12|, the GW energy emission 
drops from the (relatively large) level of the quadrupole inspiral by about an order of magnitude over a band of a few 
hundred Hz. By not rolling off more gradually, th e toy waves are likely too large as the merger frequency A/kHz is 



approached. Second, the infinitesimal bandwidth (3.15) we use for the narrow-band detectors does not describe the 
configurations that will be used in realistic second generation interferometers — the actual detectors will have some 
finite bandwidth of tens to hundreds of Hz. 

The toy model captures enough of the "feel" of the real problem that we are confident this analysis should be a 
useful demonstration of how network tuning will work with realistic waves and detectors. Nonetheless, because of the 
details that have been thrown away, we emphasize that this analysis must be followed up by calculations that are not 
oversimplified. The interplay between these neglected features is likely to be rather important — for instance, the 
finite bandwidth of the real narrow-band detectors may help to detect the relatively slow shut off that may describe 
merger waves in nature. This would also be an opportunity to test out ideas such as dynamical tuning of narrow-band 
detectors. 

IV. RESULTS 

We now test the network tuning algorithm on the toy model. We implement tuning by numerically finding the 
minima of the anticipated variance a^^^^ [which is itself built from a numerical construction of the antipicated dis- 
tribution, using Eq. ( 3.1§| ) in Eq. ( |2.9| )]. This was done with Powell's multidimensional line minimization technique. 
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using Brent's method for the successive hne minimizations (cf. Ref. This approach works very well at finding 

local minima. If the number of minima is not too large, one can simply compute the set of all local minima, then find 
the global minimum of cr^nt in that set. This is adequate for the toy model provided that the number of narrow-band 
detectors is not too large. It may not work nearly so well when more realistic waveforms and noise models are used; 
robust global minimization methods (e.g., simulated annealing) might work better in this case. 

Before exami ning merger measurements in detail, we discuss at some length the detector networks that tuning 
produces (Sec. iVA). The toy model is simple enough that we can easily understand why this algorithm tunes the 
network to the parameters that it chooses; we expect that the intuition developed here should apply more generally. 
We then look at merger measurement in Sec. [VB. We begin with a uniform prior distribution for A (imagining that, 
aside from an upper and a lower bound, we are completely ignorant of A's value) and show that repeated measurements 
of NS-NS merger drive the distribution to gradually peak on the true value A. We look at this process at two SNR 
levels (pinsp = 10 and 30), and consider measurement using 1 and 2 narrow-band detectors. 



A. Optimal detector networks 

Since our network tuning procedure works by minimizing danti it is useful to understand this quantity's properties. 
Figure ^ shows the anticipated variance when measurements are made with a single narrow-band detector. We assume 
a uniform prior distribution for A between A = 0.5 and A = 1.5, and show results for inspiral SNR pinsp = 10 (upper 
panel) and 90 (lower panel). 

In both cases, we see two minima: one is in the vicinity of 600 — 800 Hz, the other is near 1200 — 1400 Hz. For 
Pinsp — 10, the lower frequency minimum is the global minimum; for pinsp = 90, the reverse holds. This interesting 
result has a simple explanation. As we move to higher frequencies, the broad-band detector's shot noise grows 
increasingly more important, with a spectrum proportional to This makes it advantageous to place a narrow-band 
detector at high freque ncy, where it can compensate for the broad-band detector's degrading performance. However, 



our assumed waveform (3^) gets very weak at high frequencies: the amplitude (prior to the cutoff enforced by the step 
function) decays wi th fre quency as f~'^^^- The SNR measured by a high-frequency narrow-band detector is potentially 



very small [cf. Eq. ( 3.22 )], and we are more likely to get a "false positive" — an incorrect signal measurement due to 
noise at high frequencies. These two phenomena compete. Which of these two dominates in the end depends upon 
the overall signal amplitude (here set by the inspiral SNR). 

At low SNR, the need to avoid false positives tends to win out — the low frequency minimum at / ~ 630 Hz is the 
global minimum of cr^^t. As the signal strength is increased, the danger of false positives decreases. At pinsp — 90, 
the tuning algorithm places the narrow-band detector at / ~ 1380 Hz — the signal is strong enough that the tuning 
algorithm chooses to compensate for the degraded sensitivity of the broad-band detector. 

This pattern of minima location for uniform prior distribution — a high frequency minimum near / ^ 1200 — 1400 
Hz, a low frequency minimum near / ~ 600 — 800 Hz — holds as we increase the number of narrow-band detectors. 
Table | shows the narrow-band detector frequencies that minimize af^t for inspiral SNR 10, 30, and 90. Note that 
for pinsp = 10, the tuning algorithm always places the narrow-band detectors at low frequency, even for Nn — 5. At 
higher SNR, it tends to distribute the detectors more evenly. 

It's worth noting at this point that we assume the noise statistics of all our detectors is Gaussian, and hence our 
discussion of the likelihood of "false positives" is based on what may be a somewhat optimistic assumption of how the 
detectors will behave. (Since we assume that an inspiral signal has been detected and that the merger waves rapidly 
follow, Gaussianity may not be terribly optimistic — non-Gaussian events which violate this assumption will have to 
occur in a very narrow, specific time window.) Using multiple detectors in a similar band has the additional benefit 
of serving as checks on one another; it is difficult to know how one's instruments will behave until they have been 
studied and their statistical behavior is known. 



B. Measurement sequences 

Next we examine measurement sequences, following the algorithm outlined at the end of Sec. ||. As already 
mentioned, the key element of this algorithm is updating the priors after each measurement: we replace the prior 
distribution po{X) with the posterior distribution built from the measured data, P(A|A). In this simulated sequence 
of measurements, we assume for simplicity that each measurement in the sequence has identical SNR; we also (as 
discussed in Sec. ||) use a posterior distribution that has been ensemble averaged with respect to the noise. A more 
realistic calculation would randomly pick a noise realization from the distribution of possible noise functions, and 
would allow the SNR to vary from measurement to measurement in accordance with estimates of the distribution of 
binary inspiral events. 
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We first look at how well detector networks can measure A for moderate inspiral SNR (pinsp = 10) and with a 
single narrow-band detector. We show a measurement sequence for A = 0.8 in Fig. ^. The upper left hand panel 
shows the posterior probability that would be obtained if no narrow-band detectors were used for the measurement. 
It is practically flat — there is essentially no information about A in this distribution. The other panels show 
how the distribution evolves as multiple measurements are made using a narrow-band detector. The distribution 
gradually becomes peaked as more and more measurements are made. Network tuning "zooms-in" on the true merger 
parameterization A, placing the narrow-banded detector where it most effectively measures the merger waves. After 
9 measurements, the distribution is clearly peaked around A; it is markedly peaked after 15 measurements. Such a 
strongly peaked distribution indicates that A has been measured very accurately — the narrowness of the distribution 
means that the error in A's determination is very small. 

When A = 1.2, it takes many more measurements before the distributions become peaked; this is shown in Fig. |^. 
This is not surprising — the merger features are at 1200 Hz in this case, not 800 Hz, and so the expected narrow-band 
SNR is smaller by a factor of (800/1200)''/^ = 0.6. Because of the low SNR, the tuning algorithm tends to avoid 
placing the detector at high frequencies. The algorithm walks the narrow-band detector up to high frequency fairly 
slowly, doing so at a rate that ensures we have ruled out a low value of A before "testing the waters" at high values. 
After 25 measurements, we have narrowed down to a relatively broad distribution centered at A ^ 1.1; after a total 
of 50 measurements, the distribution is strongly peaked near A. 

Figure ^ shows a sequence equivalent to that used for Fig. ^, but now using two narrow-band detectors. The 
convergence onto the true parameterization is quite a bit quicker in this case. With two narrow-band detectors, the 
tuning algorithm has much more freedom to choose how it zooms in on the merger waves than when only one is used 
— the detectors can be placed at both ends of the merger frequency space and slide inward, or both can be placed at 
one end. In the case shown, the distribution found after 28 measurements with two narrow-band detectors is nearly 
the same as it was after 50 measurements with one narrow-band detector. 

The posterior distribution becomes peaked near A much more quickly if the measurements include several events at 
higher SNR. Figure ^ shows the convergence of the distribution for A when each measurement is at pinsp = 30 rather 
than 10. There is much less danger of a false positive at this SNR level, so the algorithm is able to zoom in on the 
true parameterization much more aggressively. As a consequence, the distribution is sharply peaked after only nine 
measurements — A is measured quite accurately relatively quickly at this SNR. The convergence is even quicker when 
A — 0.8: although we do not show the corresponding plots, we find that the posterior distributions become quite 
peaked about A after only 4 measurements with pinsp 30. 



V. CONCLUSIONS AND DIRECTIONS FOR FURTHER WORK 



The results of the previous section show that network tuning works effectively, at least within the context of the 
toy model: in all cases, we find that after some number of measurements the probability distribution describing what 
is known about A will become sharply peaked about A. The number of measurements that are needed depends quite 
a bit on the particulars of the waveform and of the detector network used. In particular, we can draw two conclusions 
that are not at all surprising: 

• Multiple narrow-band detectors helps the measurement process greatly. In particular, the improvement in going 
from one narrow-band detector to two detectors can be significant. 

• The rate at which measurements converge onto an accurate description of the merger waves depends strongly 
on those waves' measured SNR. These can be seen in the drastic difference in the convergence of measurement 
sequences for A = 0.8 and A = 1.2 at fixed inspiral SNR (Fig. ^ versus Fig. ^ and in the drastic difference seen 
when the wave amplitude is increased by a factor of 3 (Fig.^ versus Fig. 

From the first of these conclusions, we advocate investigating the possibility of running at least two detectors in 
the world-wide GW detector network in a narrow-band configuration. For example, LIGO already has 3 running 
interferometers (a 4 km broad-band detector and a 2 km broad-band detector at Hanford, Washington, plus a 4 km 
broad-band detector at Livingston, Louisiana). Room has been made in the facilities for additional interferometers, 
though cost is likely to limit the number that can actually be installed. If a total of 4 interferometers can be used in 
the LIGO facilities, it may be worthwhile to reconfigure one of the broad-band detectors as a narrow-band instrument 
in order to search for merger waves, assuming the loss of a broad-band detector would not seriously impact other 
science goals (e.g., if other detectors worldwide are able to fill the gap). 

Second, we strongly advocate continued theoretical efforts to model and understand the properties of merger waves. 
Of particular interest are robust characteristics such as spectral breaks and features that should be measureable 
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without needing detailed models of the waves' phasing. This kind of understanding will make it possible to choose 
a parameterization of the waves, similar to our parameter A, that leaves a strong mark on the waveform and can 
be measured reasonably well. A robust theoretical foundation for the merger waves will make it possible to choose 
our priors and configure our network usefully so that measurements will teach us about the merger waves relatively 
quickly. 

As we have repeatedly emphasized, this analysis should be considered a first, proof-of-concept presentation of how 
network tuning can work. We believe it is very important that, having presented the principles and a simple example 
of how they work, this effort be followed by a detailed foUowup analysis that uses a description of the waves and the 
detectors that is more sophisticated than our toy model. In wave modeling, we must investigate merger waveforms 
that are not as trivially simple as those used in the toy model. In principle, this is not too difficult even now — some 
groups have already produced examples of merger waves that show the influence of the dense matter EOS |ll|, p^ . 
Those waveforms could be built into this analysis without too much difhculty 

In the detector description, it is very important that the true bandwidth behavior of the narrow-band interferom- 
eters be included — our infinitesimal bandwidth description is clearly inadequate. Real detectors will be far more 
complicated than the toy description given here, and many of these complications are actually coupled (for example, 
changing the center frequency of a narrow-band detector impacts its bandwidth) . A realistic assessment of how well 
network tuning can work and how it should be implemented must take into account these various coupled complica- 
tions. It would also be useful to include acoustic narrow-band detectors in this analysis pl[ , in order to assess what 
role they could play in concert with the broad-band interferometric detectors. 

It would not be too surprising if including a non-zero bandwidth detector in this description improved the perfor- 
mance of network tuning — finite bandwidth detectors will sample a moderate span of frequency, and can thus look 
for merger power more broadly than do the toy detectors considered here. This comes at a bit of cost: the broader 
the bandwidth of the detector, the less amplitude sensitivity it has at its center frequency. A combination of very 
narrow-band and moderately narrow-band detectors may turn out to provide the best of both worlds. It would also 
be interesting to explore how well dynamical tuning might work: if it is possible to detect an inspiral in real time 
(plausible for strong sources measured in advanced detectors), we may be able to adjust the center frequency of the 
narrow-band interferometer to follow the gravitational waves as they evolve through the late inspiral and merger. 
This could significantly speed up convergence of the parameter distributions. 

Other simplifications that we have used here should also be removed in a follow-up analysis. In particular, in 
simulating a sequence of measurements we have assumed that all measurements are at some fixed SNR. This is 
obviously incorrect. In setting the measured SNR, it would be more appropriate to assume a uniform distribution in 
volume out to some distance (say 500 Mpc), and to set the distance to each source according to that distribution. 
Then, each measured event should have its angular function [defined in the text following Eq. ( 3.16| )] taken from 
the appropriate distribution (cf. Ref. jl^). Also, when we construct the posterior distribution following a simulated 
measurement, it would be much more appropriate to construct a noise instance n and hence simulate the data stream 
g, rather than ensemble averaging. 

We advocate undertaking this kind of follow-up analysis soon, so that these algorithms are well-understood when 
it become possible to actually construct the relevant detector networks, and so that what is learned from them can 
impact the design of future detectors. 
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APPENDIX A: EVALUATION OF THE INTEGRAL ^?7^ 



In this appendix, we evaluate the integral (2/7), and thereby calculate the posterior probability distribution 
^post(A|A). First, consider a single detector: the posterior probability in this case is 



fpost(A|A) =/Cpo(A) / exp 



(n|n) 



exp 



1 



{Sh\6h) - {n\dh) - - {n\n) 



Vn 



(Al) 
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where the "one-detector" inner product is 



(A2) 



Following Finn we evaluate the functional integral (Al) by treating the datastream in the limit of discrete 
sampling. Take the total span of data to be of finite duration T . Let the sampling time be At, and define 



U = lAt, i = 0,...,Af , 

0,1 = a{U) , 



(A3) 
(A4) 

(A5) 

Let a be the vector whose components are a^, an d le t C„ be the matrix whose components are Cn,ij- Note that 
C„ is a real, symmetric matrix. The inner product (A2) becomes [cf. (it), Eq. (2.20)] 



The function C„(t) is the noise autocorrelation, 

C„(r) =E{n(i)n(t + r)}„ , 

related to the spectral density by 

/■OO 



{a\b) = lim a • 
At ^ 



(A6) 



The probability distribution (Al) can therefore be written 



^'post(A A) = hm pq{\) / cxp 
At ^ 

r — ^ GO 



-n • C,;i • n - n • C"! • Sh. - -Sh. ■ C"! • Sh. 



dn 



(A7) 



This integral is of a form commonly encountered in path integral quantization, and evaluates to 



-Ppost(A|A) = lim po(A)exp 
At ^ 

= Po(A)exp 



4 " 



iSh\Sh) 



(A8) 



[We have absorbed a constant proportional to 1/ V det Cn into the prior probability po(A). 
Now consider the entire detector network: 



-Ppost(A|A) =/Cpo(A) J exp 
Recall that the "network" inner product is 



exp 



-i (^Sh\Sh) - (^n\Sh) - i {n\n) 



Vfi 



{a\b) = 4Re / df 5,(/)* [Shifty b,{f) 
Jo 



(A9) 



(AlO) 



Because the matrix Sh{f) ^ is real and symmetric, its eigenvectors are orthonormal. If the matrix of eigenvectors is 
A then A^^ = A"^, and the matrix 



is diagonal. Define further 



(All) 



(A12) 
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Using Eqs. ( [aTi| ) and (^T|), the inner product becomes 



{a\b) = 4Re / aMY AUH^'^ifrT ^',^Tbn.{f) df 
Jq 

f°° hi 
= 4Re / m)df 
Jo 

= Y^iRe r~a',{fr[Z-\f)f'b',if)df 
k -^0 

^ Y(^'k\b'k)k . 



(A13) 



The third equahty foUows from the fact that Z^^(/) is diagonal. The inner product (a'j.|6'j.)fc is identical to the inner 
product ( |A2[ ) except that [Z'f^/^]'''' is used in place of l/Sh{f). Using ( [Al3| ), Eq. ( |A9| ) becomes 



Ppost(A|A) = /Cpo(A) y exp - J] Q(5/i'fe|<5/i;, 



'k)k + {n't,\6h'k)k + {n'Mk 



Vn' 



^Po(A)]J / exp 

7. "J 



-(5K\5K)k^{n'^\5K)k~{n'Mk 



(A14) 



Since the matrix A is unitary, the functional differential element Vn' — Vn. The integral on the last line of (A14) is 
identical to the integral in Eq. (Al), except for the slightly different inner product. The result is therefore 



Ppost(A|A) = /Cpo(A)]^exp 

fc 

= /Cpo(A)exp 
= /Cpo(A)exp 



-l{5K\5K)k 



\Y.^5h',\5K), 



(A15) 



This is the result claimed in Eq. (2. 
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TABLE I: Optimal detector networks in the toy description for uniform prior probability distribution. This table shows the 
central frequencies for networks with one broad-band and Nm narrow-band detectors, for several values of inspiral SNR. 
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10 100 1000 

f (Hz) 

FIG. 1: Broad-bamd noise used in this analysis (second generation LIGO sensitivity, optimized to detect binary neutron star 
inspiral), with two examples of narrow-band interferometer noise. These examples are for resonant frequencies of 500 Hz and 
1200 Hz; each assumes a bandwidth of 50 Hz. The dashed line is a representation of the wave strain for a signal whose inspiral 
SNR is 11.7 (corresponding to an optimally oriented source at 500 Mpc, or an "average" source at 200 Mpc). 
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FIG. 2: The anticipated variance for measurements with a single narrow-band detector at inspiral SNR pinsp = 10 and pinsp = 90. 
At relatively low SNR, the tuning algorithm places the narrow-band detector at / ~ 620 Hz to avoid false positives, since it 
knows the signal will be extremely weak at high frequencies. At high SNR, the algorithm places the detector at / ~ 1400 Hz 
— the risk of a false positive is much lower, so the algorithm chooses the narrow-band detector to counter the rapidly growing 
broad-band noise. 
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FIG. 3: The evolution of the posterior probability distribution, measuring coalescences with pinsp = 10, A = 0.8, and a single 
narrow-band detector. After each measurement, we update the prior distribution, PoC-^) = -PpoBt(''^l'^)' then retune the 
network (minimizing cr^nt). In this way, multiple measurements drive the network to "zoom in" on a distribution that is peaked 
about the true value A. In the case shown here, the distribution starts to become peaked after about 9 measurements, and has 
a pronounced peak after 15 measurements. 
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FIG. 4: The same as Fig. ^ but with A = 1.2. The evolution of the posterior probability is qualitatively similar to the evolution 
when A = 0.8, but is much slower. At this signal strength, the tuning algorithm tends not to sample the high-frequency end 
of its parameter space, in order to reduce the probability of false positive measurements. Because the merger occurs at high 
frequency in this case, the algorithm takes quite a while to find the merger power. 
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FIG. 5: The same as Fig. ^ but using two narrow-band detectors. Convergence onto a distribution peaked near A is still slow, 
but not as bad as when a single detector is used. 
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FIG. 6: The same as Fig. ^ but at higher SNR: pinsp = 30. SNR has a marked effect on the rate at which the distributions 
converge onto A. The algorithm does not need to worry about false positives so much in this case, and so it is quite happy to 
assign the single detector to high frequencies when necessary. As a consequence, the distribution is quite peaked after only 9 
measurements — in marked contrast to the 50 measurements needed when pinsp = 10. 



